International  Journal  of  Environment,  Agriculture  and  Biotechnology,  5(3) 
May-Jun,  2020  /  Available:  https://iieab.com/ 


Conception  and  sizing  of  an  industrial  waste 
water  treatment  station  equipped  with  biogas 
digester,  UV  reactor  and  reverse  osmosis  unit 

Elmissaoui  S1,  Sellami  MH1  2*,  Taher  O1 

'Laboratory  of  soft  material  and  electromagnetic  modelling.  Physics  Separtment,  FST,  ,  University  Tunis  El  Manar 

2High  School  of  Engineering  Medjez  El  Bab,  IRES  A,  Jendouba  University,  Tunisia 

*Corresponding  Author 

ORCID  ID  :  0000-0003-2019-4423 


Abstract —  The  reuse  of  the  treated  waste  water  in  all  fields  (agricultural,  industrial,  domestic)  is  an 
orientation  asked  by  all  intervening  parties  in  the  water  sector  in  the  purpose  to  protect  the  environment 
and  combat  the  scarcity  of  water  resources.  This  valorization  cannot  be  done  without  making  the 
appropriate  treatment  in  order  to  respect  the  norms.  In  this  work,  we  are  developing  a  general  method  to 
conceptualize  a  treatment  station  equipped  with  biogas  digester,  UV  reactor  and  Reverse  Osmosis  (RO) 
unit.  The  modeling  approach  developed  permits  to  test  and  to  ameliorate  the  efficiency  of  each  cited 
technique.  The  main  basic  equations  used  are  those  of  Maxwell,  mass  transfer,  radiation  transfer  and 
energetic  transfer.  We  have  made  an  analytical  resolution  of  those  equations  and  we  have  expressed  the 
necessaries  formulas  for  the  design  of  each  technique.  Results  from  case  studies  were  presented. 

Keywords —  Wastewater  valorization,  Biogas  digester,  UV  reactor,  Reverse  Osmosis. 


I.  INTRODUCTION 

Actually,  at  international  level,  the  scarcity  of  water 
resources  becomes  a  real  problem  threatening  human  well¬ 
being  and  humanity  developing.  That’s  why,  valorising 
the  unconventional  water  resources  (waste  water,  salt 
water,  sea  water),  in  all  sectors  of  human  activities, 
becomes  a  priority  and  an  orientation  of  many  countries. 

After  adequate  treatments  the  use  of  wastewater  is 
considered  as  an  important  additional  water  resource  which 
could  be  of  great  interest  to  several  sectors  especially  in 
countries  with  a  semiarid  climate.  Also,  the  reuse  of  waste 
water  after  treatment  can  protect  the  environment  from 
their  reject  in  their  raw  state.  However,  the  use  of  such 
water  without  an  adequate  advanced  treatment  can  affect 
human  health,  plants,  soils  conditions  and  the  underlying 
groundwater.  Introducing  the  membrane  technique  like 
reverse  osmosis  in  the  treatment  chain  and  disinfections  of 
wastewater  by  U.V  irradiation  are  considered  as  a  credible 
alternative  to  chemical  disinfections  because  of  the  absence 
of  toxic  products.  The  functioning  of  every  device  in  the 
waste  water  treatment  chain  needs  an  important  amount  of 
energy.  Reason  why  introducing,  in  the  treatment  chain,  a 
biogas  digester  producing  energy  from  the  extracted  sludge 


can  offer  the  needed  energy  for  the  running  of  the  waste 
water  treatment  station. 

In  this  work,  to  respond  to  requests  asked  by  all 
intervening  in  the  non-conventional  water  resources,  we 
are  developing  a  global  model  permitting  to 
conceptualize  a  waste  water  treatment  plant  equipped 
with  advanced  treatment  technologies  (UV  reactor. 
Biogas  digester.  Reverse  Osmosis  unit),  to  size  each 
device/unit  and  to  test  its  efficiency.  The  main  basic 
equations  used  are  those  for  mass  transfer,  radiation 
transfer,  energetic  transfer  and  Maxwell  A  case  study  for 
the  waste  water  issued  from  a  yeast  industry  installed  in 
Tunisia  will  be  detailed. 

H.  THEORETICAL  BACKGROUND 

As  we  said  above,  valorizing  the  waste  water  in  different 
sectors  for  all  human  activities,  necessitates  a  high  level 
treatment  techniques  characterized  by  an  advanced 
efficiency  rate.  To  elaborate  a  model  permitting  to  make 
an  efficient  conception  of  each  technique,  we  have  to 
consider  the  physical  phenomena  that  take  place  inside.  In 
fact,  for  the  reverse  osmosis  unit,  we  consider  the 
membrane  as  a  polarizer  that  diffracts  particles  to  eliminate 
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by  polarizing  their  electric  and  magnetic  fields.  So  the 
efficiency  and  the  size,  of  the  membrane  depend  on  its 
electromagnetic  properties,  on  the  hydraulic  flow  through 
it  and  on  the  amount  of  electric-magnetic  fields  exchanged 
(Elmissaoui  et  al  2020submitted).  The  UV  reactor,  serves 
to  deactivate  microorganisms  (microbes,  bacteria, 
vimses. . .)  by  contact  with  an  UV  radiation  field  emitted  by 
a  UV  lamp.  The  UV  radiations  spread  in  a  hemispherical 
angle  so  they  touch  all  the  permanent  flow  of  wastewater 
passing  through  the  reaction  medium.  Analysing  the 
efficiency  of  a  U.V  photon  reactor  and  sizing  all  its 
components  are  needed  to  specify  the  amount  of  radiation 
used  to  achieve  total  or  partial  destruction  of  micro¬ 
organisms  and  to  quantify  the  hydraulic  flow.  To  realise 
this  goal,  we  have  to  determine  the  UV  incident  radiation 
profiles  and  local  volumetric  rate  of  energy  absorption 
profiles  (LVREA),  both  as  a  function  of  the  spatial  position 
inside  the  reactor. 

For  the  bio-digester,  we  will  model  and  conceptualize  a 
biogas  reactor  that  permits  both  to  produce  methane  from 
the  extracted  biomass  (sludge  from  treated  wastewater)  and 
to  cogenerate  from  it  an  important  amount  of  energy  as 
electricity  and  heat.  Knowing  that  one  cubic  meter  of 
biogas  (60%  CH4,  40%  CO2)  can  give  heating  value  as  20- 
25  MJ  which  corresponds  to  about  half  a  liter  of  diesel  oil 
(with  a  calorific  value  of  about  6  kw/m3). 

We  will  apply  principally  the  Maxwell  equations,  the 
radiation  transfer  equations,  the  mass  balance  equations 
and  the  energy  balance  equations. 

III.  BASIC  EQUATIONS  FOR  THE  THREE 
COMPARTMENTS  OF  THE  MODELING 
APPROACH 

3.1  Modeling  approach  for  the  reverse  osmosis 
compartment 

Because  of  its  elctromagmetic  properties,  the  particle  in  its 
displacement  inside  the  waste  water  was  treated  as  a  wave 
that  interferes  with  the  others  waves  representing  the  rest 
of  particles  in  the  wastewater  (Elmissaoui  et  al  submitted 
2020).  The  filtering  membrane  is  considered  as  a  polarizer 
that  will  polarize  the  electromagnetic  field  of  every 
component  in  the  wastewater.  The  capacity  of  a  membrane 
to  eliminate  particles  that  the  end  user  of  the  treated 
wastewater  doesn’t  tolerate  is  directly  linked  to  the 
electromagnetic  fields  of  the  transmitted  and  reflected 
waves  of  every  particle  through  the  membrane.  By 
applying  the  Maxwell  equation  and  the  continuity  equation 
we  have  expressed  the  electric  and  magnetic  field  of  every 
particle  in  the  wastewater  (Haralambos  1998,  Hirovuki 
T.2001,  Karine  et  al  2017,  Kiselyov  et  al  2019,). 


The  incident  electric  field  (  E’%?1,  E™)  in 

respectively  the  three  directions  (ox,  oy,  oz)  and  for  each 
particles  m  (par,  m)  in  the  wastewater  can  be  expressed  as 
follow: 


ElZ?  =  AvcOs0^r’meikr’mr 


x.inc 


gpar.m  _  ^  „ifc: 


eL 

par.m 


y.inc 


-.par.m  _ 


A„sin0 


par ,m  ikPar,mr 


The  reflected  electric  field,  ( ,  ^re/™)  for 
each  particle  m  in  the  wastewater  for  the  three  direction 
(ox,  oy,  oz)  : 

j-.par.m  n  npar,m  ikpar’mr 
K.ref  =  Bvcos0(.ef  e 

Fpar,m  _  R  pikP^mr 
hy,ref  ~  Bh «  1 


<rT  =  B-Sin° 


par.m  ikPar’mr 
ref  e 


The  transmitted  electric  field,  (Eff.f1 ,  Eyfr’fL  fox 

each  particle  m  in  the  wastewater  for  the  three  direction 
(ox,  oy,  oz)  is  expressed  by  : 


ppar.m  _  r  Qpar,m  ik^ 
Bx,tra  f-vcosBtra  e 


Tjpar.m  _ 

y.tra 


=  Che 


ik% 


>"L  —  r  Qinf) 
uz,tra  ~  ^v^L,iUtra 


par ,m  ikij 


The  incident  magnetic  field  in  unity  of  magnetic  excitation 
(  HyZ‘™’  HzZ'f)  in  respectively  the  three 

directions  (ox,  oy,  oz)  and  for  each  particles  m  (par,  m)  in 
the  wastewater  can  be  expressed  as  follow: 


H 


par,m  _ 

x,inc 


=  vfpl£1Ahcosdl 


par,m  iktr 


jrpar.m  _ 

ny,inc 


-vviEiA„e 


ifcr 


H: 


par.m  „  .  npar,m  ik?ar'mr 

=  V'p&AbSinOL^  e  1 


3.2  Modeling  approach  for  the  UV  reactor 
compartment 

In  this  part,  we  have  proposed  equations  to  predict  the 
distribution  of  UV  radiation  inside  an  open  channel  photo¬ 
reactor  (used  to  purify  a  secondary  waste  water)  in  the 
absence  of  dispersion  effects  (  reflections  and  refractions) 
and  as  a  function  of  the  spatial  coordinates  of  the  system. 
By  solving  the  radiation  transfer  equation  we  have 
expressed  the  radiation  flux  density  and  the  local 
volumetric  rate  of  energy  absorption  profiles  (LVREA)  at 
any  point  inside  the  open  channel  reactor(Santarelli  1983, 
Alfano  et  al  1986,  Lopez  and  Semel.  1999,  Sellami  et  al 
2003)  .  The  system  parameters  for  the  radiation  transfer 
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equation  were  taken  from  experimental  measurements 
(Moreno  et  al  1997,  Tosa  and  Hirata  1999,  Abdennaceur  et 
al  2000,  Sellami  et  al  2003,  Sellami  et  al  2019). 

The  equation  obtained  may  be  formally  integrated  from  a 
point  of  entrance  at  the  interior  reactor  wall  (  X  =  Xk )  to 
the  point  under  consideration  (  X  =  x  )  to  give  the 
following  expression  for  the  incident  UV  radiation  at  any 
position  inside  the  reactor: 


Iv(x,  9, 0,  t)  =  I0:V(9, 0,  t)exp 


s 

J  Kv(S,t)dS 

S=SR 


Inv is  the  spectral  specific  intensity  of  radiation  having  a 
frequency  v  and  a  direction  of  propagation  Q  at  position  S 
and  time  t.  The  spectral  volumetric  absorption  terms  is  Kv  . 


In  the  most  general  case,  radiation  may  be  arriving  at  a 
point  inside  the  media  from  all  direction  in  space.  This 
radiation  must  be  absorbed  by  a  material  point  in  space 
containing  the  different  sorts  of  bacteria  (an  elementary 
reacting  volume).  Consequently,  the  amount  of  spectral 
incident  radiation  really  used  to  eliminate  the  bacteria  is 
given  by  an  integration  over  the  solid  angle  O  for  all 
possible  situations.  The  spectral  local  volumetric  rate  of 
energy  absorption  (LVREA)  is  formulated  by  an  integral  of 
all  the  radiation  attenuations  over  the  space  of  reaction  and 
it  was  given  by  (Lopez  and  Semel.  1999,  Sellami  et  al 
2003): 


r‘t>2  r@2 

e£(x,  t)  =  Kv(r,  t)  I  d<t>  I  d6  sinO  I0  v(8,<P,t) 

J  Jr, 


*  exp 


[  Kv(S,t)dS 
Js=SR 


Hence,  we  can  say  that  the  local  volumetric  rate  of  energy 
absorption  LVREA  (eva),  represents  the  amount  of  photons 
(in  units  of  energy  for  a  given  frequency  interval)  that  are 
absorbed  per  unit  time  and  unit  reaction  volume.  So  it  can 
be  defined  as  the  rate  of  radiant  energy  absorbed  by 
particles  and  it  permits  to  determine  the  absorption 
efficiency  of  our  reactor.  Hence,  knowledge  of  the  LVREA 
spatial  distribution  is  a  key  variable  for  photocatalytic 
reactor  design. 

3.3  Modeling  approach  for  the  biogas  digester 
compartment 

In  this  part  by  applying  conservation  equation  s  of  mass, 
energy  and  motion  we  have  formulated  the  pressure  at  any 
point  inside  the  open  channel  UV  reactor,  the  optimal 
hydraulic  and  biomass  retention  times,  the  amount  of 
methane  extracted  from  the  biomass  produced  and  finally 


the  sizing  parameter  o  the  reactor  (volume,  surface,  length  , 
depth,  height...) 


3.3.1  Hydraulic  modeling  or  the  biogas  digester 


The  hydraulic  pressure  inside  is  determined  by  applying 
the  Bernoulli  theorem  between  the  entering  point  and  any 
position  in  the  reactor  (Sellami  and  Marzouk  2013, 
Sellami2015).  We  can  express  it  by: 


Hi  Jsing.k 

ii  i  k 

Hi:  The  total  charge  or  pressure  at  the  knot  i  ,jun,i ■  The 
linear  lost  of  pressure  at  the  section  i ,  Lp.  The  length  of  the 
section  i ,  JSingy.  The  singular  loss  of  energy  in  the  knot  k 
The  total  charge  at  any  section  i  inside  the  installation  is 
formulated  by: 


Pi  Vf 
Hi  =  Zi+  —  +  ^~ 
pg  2  g 

ZL\  Altitude  of  the  point  I,  Pp.  Pressure  at  the  point  I,  VL : 
Speed  of  the  flux  at  the  point  I,  g:  Gravity  acceleration,  p: 
Volume  mass 

There  is  many  formulas  that  express  the  linear  and  singular 
lost  of  energy  and  they  depend  on  the  roughness  of  the 
inner  side  of  the  accessories  through  which  the  sludge 
passes  and  the  kind  of  singularity.  They  can  be  represented 
by  the  following  general  forms: 


V2 

Jsmg  ~ 


hn  = PQnDm 

Q  :  Debit  of  the  waste  water,  D:  Diameter  of  the  tube,  P: 
Coefficient  of  roughness  for  the  linear  loss  of  energy,  e: 
Coefficient  of  singularity,  V :  Speed  of  the  water  flux,  n,  m: 
Coefficients  characterizing  the  fluid  and  the  installation 


Then  we  can  formulate  the  hydraulic  retention  time  (Hrt) 
by: 


Hrt  — 


VolcUg 

Qeff 


Hrt  :  Hydraulic  retention  time,  Voliig  :Volume  of  the 
digester,  Qeff'-  Debit  of  the  effluent 

This  parameter  (HRT)  gives  an  indication  about  the  optimal 
time  during  which  we  keep  the  biomass  inside  the  reactor 
in  order  to  produce  an  important  amount  of  biogas  in  a 
short  period  (Gasparikova  et  al  2005,  Venkatesh  et  al 
2013).  We  define  also  the  biomass  retention  time  (Brt): 

_  _  T O^-dig^t-dig 

RT~  QeffAb-eff 
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Ab-dig  ■  Amount  of  the  biomass  in  the  digester  ,  Ab_eff  : 
Amount  of  the  biomass  in  the  effluent.  The  time  of  rest 
inside  the  reactor  in  order  optimize  its  efficiency  is: 


Etot-Vrod ■  Total  energy  produced,  Ameth_vrod  :  Amount 
of  methane  produced,  Dmeth  :Methane  density, 
Eprod/gr-meth'-  Energy  produced  per  gram  of  methane 


tync-f 


igrMndcDQ 
Ff  Sp  +  MndCDO 


r dec 


Generally  we  consider  that  95%  of  the  energy  restituted 
from  the  biogas  could  be  valorized  then  we  can  write: 

Eval  ~  0,95  Etot-prod 


trest :  Time  of  rest  (day),  xgr:  Rate  of  growt  (g/gday), 
Mnd-CDO  '■  Mass  of  non-degradablechemical  demand  in 
oxygen  (g /m3)Fdec  :Decadence  factor  (g/gday).  F1/2sp : 
Factor  of  half  the  speed  (mg  / 1) 

The  biogas  produced  is  collected  at  the  superior  part  of  the 
digester  and  the  optimization  of  the  amount  collected 
depends  on  three  principal  factors:  the  composition  of  the 
substrate,  the  time  of  rest  and  the  temperature  inside  the 
digester.  To  increase  the  efficiency  of  the  reactor  we  must 
propose  the  size  that  permits  to  have  equality  between  the 
hydraulic  retention  time  and  the  biomass  retention  time 
(Azimi  and  Zamanzadeh  2004,  Banu  et  al  2007). 


3.3.2  .  Energetic  modeling  for  the  biogas  digester: 


To  model  the  energetic  process  inside  a  digester  we  must 
consider  the  amount  of  methane  to  have  inside  the  biogas 
disengaged,  That  amount  of  biogas  released  is  directly 
related  to  the  quantity  of  biomass  to  restrict  from  the  waste 
water.  The  amount  of  biomass  produced  is  expressed  by 
the  following  equation  (Sellami  and  Marzouk  2013,  Lise  et 
al  2008,  Mogens  et  al  2008,  Arsova  2010): 


Epro  —  1,17  btrestE  deb  Edec 
+  QeffiMsus 


Ybio(Mn 


Mi-CDO ) 


1  +  TrK,tFd 


■  Koi) 


Evai :  V alorized  energy 

From  the  amount  of  energy  we  can  valorize  we  must 
consider  the  needed  energy  for  the  boiler  to  heat  the  mud 
which  is  generally  evaluated  at  about  30%  from  that 
valorized  and  the  pumping  energy.  We  noticed  them  by 
auxiliary  energy  and  we  can  write: 

p 

Uaux 

=  0,3  ^ 

Yb 

PPlQeff(£iJ sing.i  2!  j  Jlin,jLi-i+ 1  +  Zi^M+l) 

YtYm 

Eaux '■  Auxiliary  energy,  Yb  :  yield  of  the  boiler,  Yt:  Yield 
of  the  turbine,  Ym:  Yield  of  the  motor,/Sjn  j :  Singular  loss 
of  energy  in  the  knot  I,  Jun ,j  :Linear  loss  of  energy  for  the 
section  j  between  two  successive  knots,  Lj_i+ p.  Distance 
between  two  successive  knots,  Azt  i+1  :  Difference  of 
altitude  between  the  knots  i  and  i  +  1 

The  loss  of  energy  through  the  repartitions  of  the  digester 
can  be  by  convection  and  by  conduction.  Then  the  total 
loss  of  energy  can  be  formulated  as  follow: 

E lost,  rep  —  t  SkUk  A Tk  +  ^  VolumiCi  AT)] 
k  l 


Qeff  :Debit  of  the  effluent  —,Ybio  :  Biomass 

gr 

yiel — ,Mi_CDC0:  Initial  mass  of  the  chemical  demand  of 

gr 

oxygen  — ,  MndCD0:Mass  of  the  non-degradable  chemical 

or 

demand  of  oxygen,— ,  Bdeb  :  Biomass  fraction  as  debris 

(gr\  gr 

— J,MSUs:  Mass  of  suspended  matter  — , Mvol\  Mass  of 

volatile  suspended  matter  gr/m3 

The  amount  of  methane  produced  is  estimated  by: 

A-meth-prod  ~  Frneth l-CDO  ~  MndCDo)Qef f 

Rmeth:  Rate  of  methane  produced 

The  amount  of  energy  produced  is  formulated  by: 

Ftot-prod  ~  A meth-prodFmethFprod/gr-meth 


Fiost.rep  : Total  loss  of  energy  through  the  repartitions  of 
the  digester,  t :  Time 

Sk  :  Surface  of  the  repartition  k,  Uk  Coefficient  of 
conductivity  of  the  repartition  k  ,  ,  A Tk  :  Difference  of 
temperature  between  the  repartition  k  and  the  exterior, 
Volum.i  :  Limited  Volume  inside  the  digester,  C,  :  Thermal 
volumetric  capacity  of  the  fluid  in  the  limited  volume 
inside  the  digester,  A 7):  Difference  of  temperature  between 
the  limited  volume  and  the  exterior 

The  energetic  balance  for  the  digester  can  then  be  written 
as  follow: 

Free  ~  Fvai  —  Eaux  —  Elost  rev 
Erec:  Recovered  energy,  Eval :  Valorized  energy,  Eaux: 
Auxiliary  energy,  Elost  rep :  Lost  of  energy  by  repartition 
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3.3.3.  Modeling  for  the  assessment  of  the  reactor  size’s 

Generally  to  define  the  optimal  size  of  a  biogas  reactor  we 
need  to  know  the  quantity  of  biomass  to  attract  from  the 
waste  water,  the  amount  of  biogas  to  produce  and  our  real 
need  from  energy  (heat  or  electric  form).  After,  as  many 
authors  (Lettinga  and  Hulshoff  1991,  Lew  et  al  2004),  we 
will  use  simple  formulas  linking  the  input  parameters 
characterizing  the  waste  water  flux  and  the  sizes  of  the 
reactor  to  make  the  optimal  conceptualization. 


The  total  volume  of  the  liquid  for  the  acceptable  organic 
charge  inside  the  reactor  is  expressed  by  (Duncan  .M  2004, 
Mogens  et  al  2008): 


V  Olliq-tot 


1  Mi -CDoQeff 

F  T 

Lorg 


Volnq—tot'  Total  volume  of  the  liquid  in  the  reactor  (m3), 
rorg  :  Rate  of  organic  charge  kg  /m3),  Fe  :Factor  of 
efficiency 

For  the  transversal  surface  of  the  reactor  it  can  be 
estimated  from  the  following  usual  formula: 

„  Qeff 


Jtrans-reac 


v 


sp.asc 


Strans-reac  '■  Transversal  surface  of  the 
reactor(m2),  Vsp  asc:  Velocity  of  the  ascendant  flux  ( m/h ) 

The  height  of  the  liquid  in  the  digester  is  calculated: 


„  _  V0Iu,_tot  ^  Volgaz-pro 

rlliq-dig  ~  ^  ^ 

° trans-reac  ^ dir 

Volgaz-pro  :  Volume  of  gaz  produced,  Sdir  :  Direct 
surface  of  the  section 

IV.  RESULTS  AND  DISCUSSION 
In  an  advanced  waste  water  treatment  station,  the  waste 
waters  have  to  pass  respectively  through  many  preliminary 
treatment  techniques  placed  in  chain  representing  the 
pretreatment  phase.  After,  when  the  treated  waste  water 
fulfills  a  defined  conditions  of  turbidity  it  passes  through 
the  reverse  osmosis  treatment  and  then  through  the  UV 
reactor.  The  recovered  mud  will  enter  the  bio  digester  in 
order  to  produce  the  methane  which  will  be  transformed  to 
energy  (heat  or  electricity).  That  energy  can  be  valorized 
for  the  running  of  the  treatment  station.  We  will  present  in 
this  section  results  from  applying  the  equations  and 
formulas  developed  above  for  every  compartment  alone. 

4.1  For  the  reverse  osmosis  compartment 
The  final  equations  were  analyzed  and  resolved 
analytically.  Firstly,  we  have  estimated  the  evolution  of  the 
electric  field  and  magnetic  field  as  function  of  the  position 
through  the  membrane  for  the  sodium  ion. 

Secondly,  we  have  estimated  the  evolution  of  the  electric 
field  for  different  incidence  angles  and  as  function  of  the 
positions  through  the  membrane  for  the  sodium  ion.. 
Figure  n°  1  presents  the  evolution  of  the  electric  field 
through  the  membrane  for  the  incident  wave  describing  to 
the  sodium  ion  motion. 
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Fig.l:  Evolution  of  the  incident  electric  field  through  the  membrane 
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The  variation  of  the  magnetic  field  as  function  of  the  depth  through  the  membrane  for  the  incident  wave  representing  the 
sodium  ion  in  its  displacement  is  in  figure  3. 
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Fig.2:  Evolution  of  the  incident  magnetic  field  through  the  membrane 


For  both,  the  electric  and  magnetic  incident  fields,  we  can 
see  clearly  that  the  global  trends  depend  on  the  position 
inside  the  membrane.  Also  amplitudes  and  distances 
between  successive  extremes  vary  through  the  membrane. 
These  findings  clearly  show  the  influence  of  the  membrane 
structure  and  texture.  Then  we  can  consider  the  modeling 
approach  as  a  tool  to  diagnose  the  membrane  state  and  to 
test  its  efficiency.  In  fact,  we  have  considered  in  this  work 
every  particle  in  the  wastewater  in  it  displacement  as  a 
wave  that  interacts  and  interferes  with  others  waves 


(particles)  and  with  filtering  membrane.  The  most  efficient 
membrane  is  the  one  through  which  the  interference 
intensity  is  minimal.  The  modeling  approach  proposed  here 
permits  to  evaluate  this  phenomenon  by  considering  the 
membrane  structure  and  texture. 

4.2.  For  the  UV  reactor  compartment 
Figure  3  and  figure  4  present,  respectively,  the  resulting 
radial  profiles  and  the  local  volumetric  rate  of  energy 
absorbed  (LVREA)  inside  an  open  channel  photo  reactor 
containing  wastewater 


Fig.  3:  Radial  profiles  of  the  incident  UV  radiation  energy  density  (Theoretical  values  ■  Experimental  values  ) 
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Fig. 4:  Radial  profiles  of  the  LVREA  (Theoretical  values 


■  ,  Experimental  values  I  ) 


For  figure  3,  we  can  say  that  increasing  the  radial  distance, 
the  value  of  (q>)  decreases  because  we  are  moving  away 
from  the  radiation  source.  However,  it  is  interesting  to 
notice  that  the  curve  goes  through  a  maximum.  This  is  the 
result  of  in-scattering  contributions  in  the  region  closer  to 
the  inner  reactor  wall  where  intensities  have  their 
maximum  values.  From  that  point  on,  geometric  effects 
supersede  the  influence  of  in-scattering,  and  the  incident 
radiation  decreases  thereafter.  When  the  particle 
concentration  (bacteria)  is  increased,  this  phenomenon 
becomes  more  important.  At  both  ends  of  the  reactor,  this 
profile  shows  a  continuous  and  important  decrease.  This 
result  is  not  related  to  the  reactor  performance;  it  is  just  the 
direct  effect  of  the  boundary  condition  (Sellami  et  al  2003, 
Sellami  et  al  2019). 

For  figure  4,  we  can  see  that  the  absorbed  energy  is 
maximal  for  the  highest  bacteria  concentration  and 
absorption  is  almost  complete  at  the  end  of  the  reactor. 
Radial  profiles  indicate  a  strong  absorption  for  low  values 
of  r.  However,  the  maximum  absorption  is  not  observed  at 
r  =  0.  This  behaviour  should  be  attributed  to  the  presence 
of  a  concentrated  boundary  condition  in  a  narrow  field  of 


directions.  In  fact,  at  the  reactor  top  position,  this  profile 
seems  to  represent  just  an  exponential  decay.  When 
changing  the  location  of  the  radial  profile,  the  main 
difference  in  the  reactor  is  that  contributions  from  the  lamp 
volume  of  irradiation  become  increasingly  smaller. 

4.3.  For  the  biogas  digester  compartment 

Application  was  done  in  order  to  size  a  treatment  station 
for  a  yeast  industry  installed  in  Tunisia  (Sellami  and 
Marzouk  2013). The  debit  of  waste  water  issued  is  about 
1500  m3/day  with  a  charge  of  about  10000  mg/1  in  unity  of 
chemical  demand  in  oxygen  (CDO),  8871  mg/1  in  unity 
biochemical  demand  in  oxygen  in  five  days  (BD05),  886 
mg/1  in  unity  of  suspended  matter  and  486  mg/1  in  unity  of 
nitrogen.  Annually  we  can  reach  a  charge  of  about 
5612970  kg  for  the  chemical  demand  in  oxygen  and  a 
charge  of  about  4856872,5  kg  for  the  biochemical  demand 
of  oxygen. 

We  have  calculated  the  hydraulic  characteristics  of  the 
biogas  digester,  its  physical  parameters,  the  energy 
produced  and  its  thermo-economic  function.  The  following 
tables  recapitulate  our  results. 
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Table  .1:  Hydraulic  characteristics,  physical  parameters  and  biogas  produced 


debit  of  the  effluent  Q,  nr3/h 

1500 

velocity  of  the  ascendant  flux  V,  m/h 

0,55 

surface  of  the  transversal  section  of  the  digester  A,  nr 

113,13 

Length  of  the  surface  ,  m 

13 

Width  of  the  surface,  m 

9 

Total  volume  of  the  liquid  in  the  reactor  Vl,  m3 

905 

Height  of  the  reactor  Hl,  m 

8 

Height  for  the  collect  and  the  storage  of  the  biogas  Ho,  m 

2,5 

Total  height  of  the  reactor  Ht,  m 

10,5 

Hydraulic  retention  time  HRT,  h 

14,5 

Time  of  rest  TR,  Day 

4 

MES  concentration  in  the  biomasse  zone  Xmes,  Kg/m3 

2,7 

Amount  of  sludge  produced  Px.mes,  Kg/day 

3  082 

Amount  of  methane  produced,  m3/day 

3  877 

Energyproduced,  KJ/day 

123,27. 106 

Amount  of  CaCCL ,  Kg/day 

1  950 

Table.  2:  Results  from  thermo  economic  analysis 


Valorisable  energyGWh/year 

13 

Energy  for  digester  heating,  GWh/year 

4,34 

Energy  for  pumping  GWh/year 

0,02 

Energy  disponible,  GWh/year 

8,64 

Electrical  energy  produced  GWh/year 

2,6 

Designation 

Cost  DN 

Installation  (material,  work  force . . . 

30413 

Equipment 

87660 

Total 

118073 

Diverse  unexpected  (15%) 

17710 

Somme 

133494 

Yearly  valorized  energy 

13  GWhyearly 

Yearly  gain  of  energy  after  considering  the  energy  consumption 

8.64  GWhyearly 

Yearly  Electrical  energy  produced 

2,6  GWhyearly 

Power  of  the  motor  (KW) 

356 

Cost  of  the  motor  (DT) 

1254000 

Cost  of  functioning  DT/an 

364  000  DT/an 

Annual  gain  (DT/year) 

119  600 

Compensation  time  without  subvention 

10  ans  et  6  mois. 

Compensation  time  with  a  subvention  of  30  to  40% 

6  ans  et  4  mois 
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The  system  we  are  conceptualizing  in  order  to  cogenerate 
electricity  and  heat  from  the  biogas  produced  after  treating 
the  waste  water  issued  from  yeast  industry  is  formed  from 
three  important  parts:  biogas  collector  (gas-liquid-solid 
separator,  dropping  network),  tank  of  sludge,  alimentation 
system. 

The  thermo  economic  evaluation  (Erlach  et  al  1999)  of  the 
biogas  produced  from  the  two  stations  shows  that  we  can 
realize  a  gain  of  about  119600  DT/year  for  the  yeast 
industry  with  a  compensation  time  of  the  cost  of  reactor 
installation  between  6  and  10  years. 

V.  CONCLUSION 

We  have  detailed  in  this  work  a  modeling  approach  that 
permits  to  conceptualize  a  wastewater  treatment  plant  to 
valorize  the  wastewater.  The  plant  is  equipped  with  three 
advanced  treatment  techniques:  a  reverse  osmosis  unit,  a 
UV  reactor  and  a  biogas  digester.  We  have  presented  the 
necessary  equations  to  conceptualize  every  device  and  to 
size  it.  For  every  part  of  the  model  physics  analysis  was 
detailed  and  results  from  theoretical  calculus  were  given. 
An  application  to  valorize  the  biomass  extracted  from  the 
wastewater  issued  from  a  yeast  industry  was  done.  This 
valorization  consists  in  producing  energy  from  the  methane 
recovered  and  in  using  that  energy  in  the  running  of  all  the 
plant.  A  thermo-economic  analysis  was  done.  We  can  say 
that  we  can  realize  a  gain  of  about  119600  DT/year  for  the 
yeast  industry  from  the  energy  produced  with  a 
compensation  time  of  the  cost  of  reactor  installation 
between  6  and  10  years. 
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